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Abstract 

Let X = (Xi,X2,X3) be a Gaussian random vector such that X ~ A/JO, E). We consider 
the problem of determining the matrix E, up to permutation, based on the knowledge of the 
distribution of A' m in := min(Ai, X2, A3). Particularly, we establish a connection between this 
identification problem and a geometric identification problem in the context of the theory of the 
circular radon transform. 
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1 Introduction 

A Gaussian random vector is a vector valued random variable X = (X \,..., X n ) whose components 
Ah,..., X n are jointly Gaussian. The mean vector /z, of a Gaussian random vector X = (X\ ...., X n ) 
is defined by p, = (/zi,..., p n ) T , where pi = -E(AJ) for i = 1,..., n. Moreover, the correlation matrix 
E of X is an n x n matrix defined by E := (E^), where E,y = Cov(Xj. Xj) for 1 < i,j < n. The 
statement “X is a Gaussian random vector with mean vector p and correlation matrix E” will be 
compactly written as X ~ A f{p, E). 

Suppose that A' = (XJ,..., X n ) ~ Af(0, E) and let A m j n be the random variable defined by 
A m i n = min(Ai,... ,X n ). It is clear that the knowledge of E allow us to determine the distribution 
of A m j n . On the other hand, a natural question arises: Does the distribution of A m j n determine the 
matrix E? In other words, is it possible to recover the matrix E from the distribution of A m j n ? The 
purpose of this paper is to solve this problem affirmatively, under the assumption that l'E -1 > 0, 
which in particular covers the case where the correlations are negative mm- 

This problem seems to have originated from an econometrics supply-demand problem posed to 
Anderson and Ghurye J>]. Additionally, a similar set of problems were previously known in the 
context of competing and complimentary risks [RIolITTI. These kind of problems of identification have 
been the subject of interest to many authors, including mffiMg. Particularly, in the same setting 
established above, the authors in [J] studied this recovery problem under the hypothesis pijCTi < ctj. 
In [SIITT] it was studied the case of common correlations. And, in EE EH it was studied the case 
of non-negative correlations. The novelty of our approach consists in tackling the problem via a 
geometric approach, reducing it to a recovering problem in the context of circular radon transform. 
This allows us to significantly extend previous results. Moreover, from this approach much intuition 
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is gained regarding the ‘backstage’ geometric difficulty implicit in this and similar recovery problems, 
particularly improving previous results established in [9] and m- 

This paper is organized as follows: Section [2] is devoted to develop some preliminary geoemetric 
notions. Section [3] contains the proof of the main result of the article, while Section U presents in 
detail the proof of Lemma 0] which is a key component in the proof of the main theorem. 

2 Preliminaries 

2.1 Generalized square roots 

Definition 1 We say that N is a generalized square root of a matrix Ad (or, for the purpose of this 
article, just a square root of Ad) if NN* = M. 

Now, if M is an n x n positive definite symmetric matrix, it has an eigendecomposition PDP 1 * , 
where P is a unitary complex matrix and D is a real diagonal matrix whose main diagonal contains 
the corresponding positive eigenvalues. Thus, the matrix TV = PD 2 is a generalized square root of 
M , (where is a square root of D in the ‘usual’ sense). 

Note that, given an orthogonal matrix O and any square root TV of Ad, the matrix NO is also a 
square root of M. Moreover, each square root for Ad can be obtained by right multiplication of N with 
an orthogonal matrix. In other words, the orthogonal group acts transitively by right multiplication 
on the set of square roots of Ad. 

Notice that a square root N of a positive definite matrix Ad is a natural change of variables between 
the Hilbert space induced by the internal product (■ u,v) M := uf Adv and the Euclidean Hilbert space 
endowed with the usual internal product (u,v) := idu, in the sense that (u,v) M = ( Nu,Nv ). A 
particular use of this fact appears when we consider a Gaussian random vector X ~ A/"(0, E). In this 
case we have that X = NU, where N is any square root of E and U ~ jV(0, I). 

We will impose an additional condition on the square roots of E: If a square root TV of E satisfies 
that TVei is a positive multiple of the vector 1 = (1,..., 1)*, then TV will be called a standard square 
root of E. Notice that, if this is the case, TVei = /c _1 l, with k := V 1*E _1 1. 

2.2 Cones and sections 

Let A be an n x n non-singular matrix. The positive cone associated with A is defined as the set 
C A ~ {u G R n : Au > 0}Q A different representation for Ca is given by 

Ca = < CiVi : Ci > 0 for i = 1,..., n 

U=i 

where, for i = 1, ...,ra, the vector Vi is a positive multiple of A~ 1 ei. The vectors vi,...,v„ are 
collectively called the directions of Ca- 

In the following let us assume that E _1 l > 0 and let TV be a standard square root of E. Notice 
that Cn (except for the origin) is completely contained in the half space {w £ K 3 : e\w > 0}. Indeed, 
if 0 ^ u £ Cn, then 


e\u = Al* (TV -1 )* TV _1 TVu = Al‘E -1 TVu = A(E -1 1 )‘TVm > 0. 

1 < denotes the componentwise order, that is for a given vectors v = (vi ,..., n n ) £ , w = (wi, ... ,id n ) £ E M n , v > w 

sii Vi > Wi, for each i : 1 ,..., n. 


2 




In particular, since ei £ Cn , the requirement E _1 1 > 0 allows us to choose the directions of 
Cn in a such way that they have the form Vi = (l,aj), i = 1 ,n. In such a case, the vectors cq 
correspond to the extremal points of the intersection of Cn with the plane x\ = 1. We will denote by 
Tjv the convex set generated by the points a7 £ R " _1 and we will call it the section associated with 
the standard root N of E. 


3 How to recover £ from the distribution of X m - m ? 

3.1 Reducing the problem 

Given a Gaussian random vector X = {X\, X 2 ,X 3 ) ~ 7V(0, E), we define as the tail distribution 
of W m ; n = mm(Xi, X%, X 3 ): 

m E (t) := Prob (X min > t), 

In the following, we will assume that E -1 l > 0. 

Lemma 2 Let N be a standard square root ofT,. Then E is uniquely determined (up to a permuta¬ 
tion), from the vertices of the section Tn and the parameter k := V l t E _ 1 l. 

Proof. Let Wj = (1,0* ,Pi),i = 1,2,3 be the vertices of the triangle Tn and let W be the ma¬ 
trix with columns uq, W 2 , 103 . Define by the relation [/zi, /j. 2 , Hz ] 1 := nW~ 1 ei and let 

A := diag(/ii,/i 2 ,^ 3 ). We assert that there exists a permutation matrix P such that 

N = P t A~ 1 W~ 1 . 

From this statement the lemma would clearly follow since we would have that 

E = P t A~ 1 W ~ 1 (W -1 )* (A - 1 ) 4 P. 

In order to prove the assertion notice that, up to ordering and scaling, the vectors w± . W ‘2 , W 3 are equal 
to the directions _/V - 1 ei,_/V - 1 e 2 , A^” 1 e 3 , thus 

(V e\,N e-2, N 63) (^cr(l)^ ? (7(l), / ^<7(2)^ J cr(2), ^cr(3)^cr(3)) 

for some permutation a and some Ai,A 2 ,A 3 > 0. This is equivalent to saying that WA'P = N -1 , 
where P is the matrix permutation associated to a and A! := diag(Ai, A 2 , A 3 ). 

Consequently, 

WIAi.Aa^s]* = WAT = WA'Pl = N^l = Kei . 

Thus [Ai, A 2 , A 3 ]* := KW _ 1 ei, implying that A = A' from where the assertion follows. ■ 


Let N be a standard root of E. Notice that for any u £ R 3 the condition Nu > 1 1 is equivalent to 
N(u — tme 1 ) > 0, and therefore {u : Nu > fl} = tne 1 + Cn- Then, we have that for t > 0, 

ms (t) = Prob (V min > t) = n {tne 1 + C N ), (1) 

where /i is the standard Gaussian measure N {0,1) in R 3 . Moreover, if stands for the Gaussian 
measure Af (0, a 2 1) , it is the case that 

toe (t) = Hi /(*«) (ei + C N ). (2) 
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Lemma 3 As t —> oo, In ms (t) ~ — ^-§~. In particular, 


2 1 ,. In ms (t) 

t = — — Iim 


2 t—t+OO 


t 2 


Proof. The rate function of the sequence of measures p a , where a — > 0, is given by I (u) 
(see nn). Therefore, 

ln/v (ei + Cn) ~ —ct ~ 2 inf I(u); 

u(z(ei -\-Cn) 

On the other hand, it is clear that 


( 3 ) 


IMI 2 A 


inf 

u£(ci+Cn) 


J( U ) = J(ei) = l 


Therefore, as t —> oo, we have that 

In (ms (f)) = ln/x 1/tK (e x + C N ) ~ 


t 2 n 2 

2 


We define the (two-dimensional) circular transform of a function / : R 2 —> K. by 


2-jt 

Rf (p) := Jf (pcosd, psinO) dO. 
o 

For a measurable set S C K 2 , we define the circular transform by 

2 tt 

Rs (p) ■= Ri(-eS) (p) = JIs(pcosd 1 psm6)de i p> 0, 
o 

where Is denotes the characteristic function of S. Notice that, since the angular measure is invariant 
under orthogonal transformations, the circular transform is also invariant under orthogonal transfor¬ 
mations. In other words if O is an orthogonal transformation, it is the case that Rf 0 o = Rf - On the 
other hand, if Rf = R g for given functions /, g , it is not necessarily the case that / o O = g for some 
orthogonal transformation (see example [5]). In spite of this, we have a positive result in this direction 
when / and g correspond to characteristic functions of triangles enclosing the origin (that is, such 
that 0 belongs to its interior). 

Lemma 4 If T is a triangle enclosing the origin , then T can be recovered (up to an orthogonal 
transformation), from Rt- 

Due to its length, we relegate the proof of Lemma |H to Section |H Although moderately technical, it 
relies on elementary geometry and basic linear algebra. The condition that the triangles enclose the 
origin is necessary, as the following example shows: 
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Example 5 Consider the triangles depicted in the figure below 


(0,4) 



(0,3) 



( 0 , 0 ) 


(4,0) 


It is clear that the triangles have the same circular transform. Therefore, for the triangles 




the circular transform is also the same. However, it is clear that they are not orthogonally equivalent. 


3.2 The main Theorem 


Now we are able to state and prove the main result of this article: 

Theorem 6 Suppose that X = [X\, X 2 , X3) ~ A/"(0, £). Let X m j n be the random variable defined 
by X m j n = min(Xi, X 2 , X3) and assume that l'S -1 > 0. Then, the distribution of X m j n uniquely 
determines £ up to permutation equivalence. More exactly, if Y = (Yi,Y 2 ,Ys) ~ 7V(0, £ 0 ) and 
Xmin = Tmin> there exists a permutation matrix P such that £q = P£P ( . 


Let us first prove the following lemma: 

Lemma 7 Let N be a standard square root of £ (see definition \7\) where 1 *£“ 1 > 0. Then Rt n is 
identifiable from ms (t) := Pr(X m ; n > t). That is, ms ( t) uniquely determines Rt n 

Proof. From eq. Q, we have that 


mu (t) = p i K te 1 + Cat) = 


(2tt) 


3/2 




e 21 


Ktei+CN 


where n 2 = 1 4 £ 1 1. Equivalently, 


+ OO 

ms (t) =- f jj e _ 2 “ 1 e _2 ^ 2+ “ 3 ^ du 2 du 3 dui, 

(27t) J J J 


fit ( ui—K,t)Tjsr 


5 










Now, a change of variables leads to 


+OO 

m s (t) = — J e ~^ x+Kt)2 Jj e~^ y2+z2) dydzdx 

0 xT n 

+OO 

Je-««e-i*’JJ e -^ dtdzdx . 


(2tt) 

v ' 0 xT n 

Therefore, if £ denotes the Laplace transform and h is the function defined by 


h(x) := e~^ 2 Jj e-H^ 2 ) dydz , 
xTn 


we have that 


+ OO 

(2t r) 3/2 e^ms (i/«) = J e~ tx e"^ Jj e~^ 2 +* 2 ) dydz 


xTn 


dx — £ (h) ( t ). 


Since the Laplace transform is injective over the functions of polynomial growth on (0, oo), in partic¬ 
ular, we have 


£ 


or equivalently, 


-l 


((27r) 3/2 e^ ms (t/nj) (v^) = h(Vte) = e x JJ e 2 ^ V+Z ^dydz, 


V2xT n 


—C 1 ^(27r) 3 ^ 2 e 2 * ms (t/ft)) (V2x) = — JJ e 2 ( V+Z ^dydz. 


\/2xTn 


On the other hand, it is clear that 

- f[ e-^ y2+z2) dydz = 2 f fe- x(y2+z2) dydz 


\plxT n 


T n 


2 JJ e xp It n (p cos 0 , p sin 0)p dO dp; 


(0,oo) X (0,27r) 
+oo 


= / e 


= £(s 0 W, 


ItAVp cose, y/p sin 9)dQ 


dp 


(4) 


where g denotes the function g{p) := J lT N {yfpcos9, y/p sind) 

o 

Therefore, it follows that 


/ 


£ 


-1 


] 27T 

(p 2 ) = J It n (p cos 9 , p sin 9) dd. 


i x 

\ V^xTn 


(5) 
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Finally, by combining Eqs. 0 and 0 we obtain that 


R Tn ( p) = C 1 1 ((27r) 3/2 e2 t2 m E (t/«)) (v^)^ (p 2 ) . ( 6 ) 

Therefore, equation 0 give us an expression for Rt n in terms of m E (t) and k. However, from lemma 
[3]we know that k is recoverable from m E (t), thus the rigth hand side of 0 is determined from ?n E (f), 
and the result follows. ■ 


Proof of Theorem 0 From lemma EJ k is recoverable from m E (t) using eq. 0. Moreover, from 
lemma [0 we can recover Rt n from to e (t) using the formula stated in eq. ©, where TV is a standard 
square root of E. Therefore, using Lemma 2] we can recover Tjv up to an orthogonal transformation, 
that is, we can recover OT/v where O is an (unknown) orthogonal transformation in R 2 . Now, we 

. Notice that OT N = T N q, so that OT/v is the section associated with the standard 

root TV := NO. Therefore, since we can recover Ty, and k, from lemma [5] we can recover E up to 
permutation equivalence. ■ 


define O = 


1 0 
0 O 


4 The circular transform of a triangle 

We can think of the circular transform as a systematic ‘scan’ that recognizes the mass at distance p 
from the origin. When the scan encounters an abrupt change of media, its smoothness is momentarily 
lost. For instance, in order to recover an acute triangle T containing the origin, by detecting changes 
in the smoothness of Rt we can determine the distance from the origin to the sides and vertices of 
the triangle and then rely on a geometric construction to recover T. However, we take a detour from 
this approach and we use a more concise tool. Namely, that the circular transform of some ‘basic 
triangles’ form a linearly independent set (lemma 0. Then, we use such basic triangles as building 
blocks to recover more complex geometric figures, in particular, any triangle containing the origin. 

4.1 Parametric form of a triangle 

Let T be a triangle with vertices x,y,z £ R 2 and such that 0 belongs to the interior of T. Also, let 
n xy ,n yz ,n xz be the closest points to the origin from the lines xtj, yt and xt respectively (we will call 
them the heights of T). We define the parametric form of the triangle T as the ordered sequence of 
distances (||a;|| , \\n xy \\ , ||y||, \\n yz \\ , ||z||, ||n X2 ||). It is an easy geometric fact that this list determines 
the triangle T up to rotation. 

As an intermediate step we require the following lemma, which states that, if the terms of the 
parametric form can be recovered by pairs, then the parametric form of the triangle can be recovered. 

Lemma 8 Consider a triangle T with parametric form (||x|| , ||n xa || , ||y|| , ||% 2 || , || 2 /|| , ||n X2: ||). Let 
ot,/3,7 , »7i, ?72,773 satisfy 


f 

(771,0), 

(772,0)4 _ 

f(IKJ,IMI), 

(IKJ, IMI), 

(IK 4 UMD 4 

1 (*72,7), 

(773,7), 

(773, a ) J 

UIK.IUMI), 

(IKzIUMI), 

(IK*II ,IMI) J 


Then, a parametric form of T is given by 

(a,r?i,/3,T72,7,%) • 
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Proof. Without loss of generality, we can assume that a = ||a;||, /? = ||t/|| and 7 = ||2||. If it is the case 
that a , ft and 7 are different numbers, then it is clear that ||n xy || = 771, \\n yz \\ = 772 and ||n X2 || = 773. 
On the other hand, if two of these numbers are equal, say a = ft 7^ 7, then necessarily ||n xy || = 771 
and {||%z||, \\n xz \\} = {772,773} where the two possible choices lead to isomorphic triangles, therefore, 
by appropriately interchanging y and z, we have that \\n yz \\ = 772 and ||n X z|| = 773. Finally, if it is 
the case that a = f 3 = 7, then we have that {Hn^H , \\n yz \\ , ||n X2 ||} = {771,772,773} and in fact any 
assignation leads to isomorphic triangles. Therefore, appropriately interchanging x , y and 2, we have 
that ||?r xy || = 771, ||n y2 || = 772, \\n xz \\ = 773. ■ 


4.2 Some properties of the circular transform 

A family of functions {f a (:r)} QgA is said to be linearly independent in an interval / if for any finite 
set F C A, YlaGF c « fa (a) = 0 for all x £ I implies that c a = 0 for all a £ F. 


For any a,b > 0 , we define the triangle 

T (a, b) := conv ^( 0 , 0 ), (a, 0 ), ^a, yjb 2 — a 2 j j 

where conv stands for the convex closure in R 2 . We also define the function <J> aj b : ( 0 , +00) —> R by 

$£ 1 , 7 ) (p) := t?T(a,b) (p) • 

Lemma 9 The set of functions { < f > a ,b} a b>0 is linearly independent in ( 0 ,+oo). 

We first prove the following lemma. 


Lemma 10 Consider an increasing sequence of real numbers 0 < ao < • • • < a m . Then, for every 
x £ (0,1/a^) and every i, j = 1,..., m, i 7^ j, it is the case that 


a i (1 - alx) 7^ a 2 ( 1 - a 2 x ) 

Proof. Let x £ (0, l/a„) be fixed. Since 0 < x < 1/a^ < 2/a 2 , for all i = 1,... ,m, we have that 
at £ ^ 0 , \j2/x^ for 7 = 0 ,... , m. The result follows from the fact that the function f(t) := t 2 (l — t 2 x) 
is increasing in (0, \j2jx'). ■ 


Proof of Lemma [9j Given 0 < ao < • • • < a m , consider the functions {tf ai (x)}"f 0 , where 

1fa (x) := (1 - a 2 x) 1/2 . 

The Wronskian of these functions in ( 0,1 /a^) is given by 


W(tjj ao ,...,ipa m ) =C 


1 

a o (! - a o a -) _1 
al m (1 - a%x)~ m 


1 

a m (! - a 2 m x Y 

(! - a? m x)~ m 


c = yr m MItt 


i =0 


(1 - a 2 x) 


- 1/2 


where 






The matrix in the expression of W (if ) ao ,..., ip am ) is a Vandermonde matrix, with generators a 2 (l — afx) 
for * = 0,..., m . It follows from LemmaQillthat these generators are all different. Therefore the Wron- 
skian W (ip ao , ■ ■ ■ a m ) is non zero in (0,1/a^J. Consequently, for any fixed interval I C (0,1/a)/), 
the functions {ip ai (*)}”!. 0 are linearly independent in I. Therefore, the functions {aiXip ai {x)}”P 0 are 
linearly independent in I C (0,1/a)/). In consequence, the change of variables x = 1/p 2 , produces 
linearly independent functions in I C (a m ,+oo), given by 


eg 

pVp 2 ~ a V 


i = 0 , 


m. 


Moreover, their antiderivatives {arccos (a^/p) + 7i}™ 0 , where 7o, ■ ■ ■ ,7m are arbitrary constants, will 
be also linearly independent in / C (a m ,+oo). 

Now, in order to prove that the functions { < h aj f,} Q b>0 are linearly independent in (0,+oo), we 
assume that 

^ ^ l » ^ 777 . 

E l=0 ^ = °^ p > °; 

where 0 < a*,o < • • • < Oi jmi for i = 0 ,..., l, 0 < bi < ■ ■ ■ < bi and all the coefficients Cij are nonzero 
(that is, we assume the existence of a nontrivial minimal dependent set). 

Now, eq. (S3 implies that (p) = 0 , whenever i < l and p £ (&z_i, 6 j), therefore 

£™/ 0 (p) = 0 , Vp G (bi-ubi) ■ (7) 

Now, since 


^a.b (p) = arccos ( a/b ) I (p < a) + [arccos (a/b) — arccos (a/p)] I (a < p < b ), 

the eq. [7] restricted to the interval (az, mi , bi) takes the form 

£”=o c Z 7 [arccos(a;j/p) - arccos (aij/bi)\ = 0, for all p G ( ai tm ,,bi ), 

which contradicts the fact that the set {arccos (aij/p) + with 7 ^ = — arccos (aij/bi), j = 

0 ,..., m is linearly independent in the interval (ai t mi ,h). ■ 


4.3 Decomposing the triangle 

In order to prove Lemma 0] we will distinguish the following three cases: 

• Case I: All the heights of T lie in the interior of the sides. That is n xy = s\X + (1 — Si)y, 
n yz = s 2 y + (1 - s 2 ) z and n xz = s 3 x + (1 - S3) z, for some si, s 2 , S3 G (0,1). 

• Case II: There is one height that does not lie in the interior of the corresponding side. Say, 
withouth loss of generality, n xy = sx + (1 — s)y, where s > 1 . 

• Case III: There is one height that lies over a vertex. Say, withouth loss of generality, n xy = x. 
We define the basic subtriangles of the triangle T as the following six right triangles: 

T\ = conv (0, x, n xy ), T 2 = conv (0, y, n xy ) , T 3 = conv (0, y, n yz ), 

T 4 = conv (0, 0 , n yz ) , T 5 = conv (0, z , n xz ) , T 6 = conv (0, x, n xz ). 

Notice that these triangles are non degenerate (that is, of positive area), except in case III, in which 
Tf is a degenerate triangle. We should also point out the following elementary facts: 

Let I a : R 2 —> {0,1} denotes the characteristic function on the set A C R 2 . Now, 


-1/2 
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• In case I, 

: T = EL-i J Ti ■ 

• In case II, 

I T = —I Tl + Ei = 2 7 T», 

and the subtriangle T\ is not congruent to any other of the subtriangles. 


• In case III, 

= Ei=2 J Ti- 


(8) 

( 9 ) 


( 10 ) 




Given a basis B = {vi}^ =1 of a vector space and an expression 

w = ^Pt j v tj ( 11 ) 

3 

where {tj}j is an arbitrary sequence in {1 ,..., d}, we say that eq. (fill) is irreducible if Et=i l a *l = 
Ej \&tj I where oq ,,ad are the coefficients of w in the basis B. 

Lemma 11 Rt belongs to the linear span of {$ a ,&} a &>0 . More exactly, if r /1 := n xy , ??2 := n yz , 
ri3 ■= n xz , 

1. In case I, 

R t = $IMI,IMI + $|M,||i/|| + $ IMI,IMI + $ IMUMI + + $ IM,||*||; 


2. In case II, 

r t = -^Umll.llxll + + $ IMI,IMI + $|MI,INb 

3. In case III, 

R T = $ MI,l|y|| + $ IMI,||y|| + ^IMI.IWI + $ IMI,IM| + 

Moreover, all of the previous linear representations are irreducible. 

Proof. By using the linearity of the circular transform, all three formulas follow directly from Eqs. 
@, © and m respectively. The irreducibility is clear in cases I and III due to the positivity of the 
coefficients. On the other hand, to show the irreducibility in case II we must show that the function 
^11*7! II >11*11 * s n °t * n the se t 

{^IlmlUMI* ^IMI,||y||’ ^||»72||,||*||i ^IMI,IM|> ^IMMMI} 1 

Indeed this is true since, in this case, the triangle T\ is not congruent to any other subtriangle. ■ 
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4.4 Recovering the triangle 

Proof of Lemma [4j From Lemma El Rt can be expressed in basis { < F a ,&} a h>0 , say 

Rt = J2 c a,b®a,b- (12) 

Moreover, depending if Yh c a,b is equal to 6,4 or 5, we can distinguish if we are in case I, II or III 
respectively. Notice also that the representation in eq. m can be expressed as: 

Rt = C$ai,bi + ^ai,& 2 + 3 >a 2 ,& 2 + < ha 2 ,b3 + 'ha 3) fc 3 + $1x3,6,. 

for some Oj, 6j, i = 1,..., 3 and where c is equal to 1, —1 or 0, depending if we are in case I, II or III. 
We claim that (b\, ai, 62,02,63,03) is a parametric representation of T. To see this, notice that in the 
cases I and II, from the uniqueness of the linear representation we have that 


f (oiA), 

(ai,b 2 ), 

( a 2 , b 2 ) , 

l I (KJ.IMI), 

(KJ.IMI), 

(IK.IUMI), 1 

1 (a- 2 ,b 3 ), 

( a 3 , b 3 ) , 

(03,^1) 

II (KJ.IMI), 

(KJ.IMI), 

(K.IUMI) J 


Then, from Lemma[U it follows that (bi,a±, b 2 ,a 2 , b 3 , <23) is a parametric representation of T. In case 
III the claim follows the same argument, by adding the missing pair (ai, aq). ■ 
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